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the Sudakov resummation, which renders perturbative expansions more reli- 
able. It is observed that perturbative QCD is applicable to A^ -^ Ac decays 



volutions of hard b quark decay amplitudes with universal Af, and Ac baryon 



,i^ . for velocity transfer greater than 1.2. Under requirement of heavy quark sym- 



metry, we predict the branching ratio B{Ab -^ Acli^) ~ 2%, and determine 



j^ . the Ab and Ac baryon wave functions. 
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I. INTRODUCTION 

Analyses of exclusive heavy hadron decays are a challenging subject because of their com- 
plicated QCD dynamics. Recently, we have proposed a rigorous theory for these processes 
based on perturbative QCD (PQCD) factorization theorems ||l],0. In this approach heavy 
hadron decay rates are expressed as convolutions of hard heavy quark decay amplitudes 
with heavy hadron wave functions. The former are calculable in perturbation theory, if pro- 
cesses involve large momentum transfer. The latter, absorbing nonperturbative dynamics 
of processes, must be obtained by means outside the PQCD regime. Since wave functions 
are universal, they can be determined once for all, and then employed to make predictions 
for other processes containing the same hadrons. With this prescription for nonperturbative 
wave functions, PQCD factorization theorems possess a predictive power. 

For semileptonic decays, the PQCD approach complements heavy quark symmetry in 
studies of heavy hadron transition form factors 0. Heavy quark symmetry determines the 
normalization of transition form factors at zero recoil of final-state heavy hadrons, up to 
power corrections in 1/M, M being the heavy quark mass, and up to perturbative correc- 
tions in the coupling constant as- While PQCD is appropriate for fast recoil, the region 
with large energy release, and gives a dependence of transition form factors on velocity 
transfer. For nonleptonic decays, PQCD is a more systematic approach compared with the 
phenomenological Bauer-Stech-Wirbel (BSW) model |@. In PQCD factorization theorems 
contributions to nonleptonic decay rates characterized by different scales are carefully ab- 
sorbed into different subprocesses, among which renormalization-group (RG) evolutions are 
constructed 0, leading to a scale and scheme independent, gauge invariant and infrared 
finite theory ^j. Not only factorizable but nonfactorizable contributions can be evaluated 
1^. The BSW model considers only factorizable contributions: two fitting parameters ai 
and 02 are associated with external and internal Pi^-emission form factors, respectively. Non- 
factorizable contributions must be included as additional parameters 0. 

The above PQCD formalism has been applied to heavy meson decays successfully. It is 



then natural to extend the formahsni to more comphcated heavy baryon decays. In [§] we 
have developed factorization theorem for the semileptonic decay A^ -^ pW, in which Sudakov 
resummation of double logarithmic corrections to the A^ baryon wave function was included, 
and a full set of diagrams for the hard b quark decay amplitudes was calculated. This is an 
analysis more complete than the work in the literature [Q. On the other hand, 6-baryons 
have been observed in experiments at LEP and at the Tevatron. Masses and decay widths 
of the lightest 6-baryons, as compared with theoretical predictions, have stimulated many 
interesting discussions and investigations |llO|- |l^] . When Run II of the Tevatron comes up 
with a vertex trigger employed, it will be expected to collect millions of 6-baryon events. 
Therefore, an intensive study of exclusive heavy baryon decays is urgent. 

Exclusive heavy baryon decays are dominated by 6 ^ c modes. In this paper we shall de- 
velop factorization theorem for the semileptonic decay A^ -^ AJu, and locate the kinematic 
region where PQCD is applicable. It will be shown that PQCD predictions for the involved 
transition form factors are reliable at fast recoil of the A^ baryon with velocity transfer 
greater than 1.2. Under requirement of heavy quark symmetry, we predict the branching 
ratio B{Ah -^ AJu) ~ 2%. We shall also determine the unknown parameters in the A;, and 
Ac baryon wave functions, which can be employed to study nonleptonic A^ baryon decays 
because of the universality. 

In Sec. II we develop factorization theorem for the semileptonic decay A5 —>■ AJu. Su- 
dakov resummation of double logarithmic corrections to the process is performed. The fac- 
torization formulas for the involved heavy baryon transition form factors and their numerical 
results are presented in Sec. Ill and in Sec. IV, respectively. Section V is the conclusion. 

II. FACTORIZATION THEOREM 

The amplitude for the semileptonic decay A^ -^ Adv is written as 

M = ^vJY{i-i5)MMp')\c^,{^-i^)b\Mp)) , (1) 



where Gp is the Fermi couphng constant, Vcb is the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element, p and p' are the A^ and A^ baryon momenta, respectively. All QCD dynamics 
is contained in the hadronic matrix element 

M,^{Ac{p')\n,{l-l,)b\Mp)), 

+Mv')[9i{q^)l'l^ - t92{q'y^''l,qu + g^iq^Hq^Mp) ■ (2) 

In the second expression A^^ has been expressed in terms of six form factors fi and Qi, where 
A(,(p) and Ac(p') are the A^ and Ac baryon spinors, respectively, and the variable q denotes 
q = p — p'. In the case of massless leptons with 

q,lYi'^-^,)ui = 0, (3) 

the form factors /s and g^ do not contribute. Since the contributions from /2 and g2 are 
small, we shall concentrate on /i and gi in the present work. 

The idea of PQCD factorization theorems is to sort out nonperturbative dynamics in- 
volved in QCD processes and factorize it into hadron wave functions. Nonperturbative dy- 
namics is reflected by infrared divergences in radiative corrections to quark-level amplitudes 
in perturbation theory. The construction of factorization theorem for the decay A^ —>■ AJu 
is basically similar to that for the decay A;, -^ pW in p|. The lowest-order diagrams for 
b ^ c decays are shown in Fig. 1, where two hard gluons attach the three incoming and 
outgoing quarks in all possible ways. We then investigate infrared divergences from radia- 
tive corrections to these diagrams. Small transverse momenta kp are associated with the 
valence quarks, such that they are off mass shell a bit. The transverse momenta kx serve 
as a factorization scale, below which dynamics is regarded as being nonperturbative, and 
absorbed into Af, and A^ baryon wave functions, and above which perturbation theory is 
reliable, and radiative corrections are absorbed into hard b —^ c decay amplitudes. 

Infrared divergences from radiative corrections are collinear, when loop momenta are 
parallel to an energetic light quark, and soft, when loop momenta are much smaller than 



the Aft baryon mass M\^. Collinear and soft enhancements may overlap to give double 
logarithms. Three-particle reducible corrections on the Af, baryon side are absorbed into 
the Afo baryon wave function. If the light valence quarks move slowly, collinear divergences 
associated with these quarks will not be pinched [^, and soft divergences are important. 
However, there is probability, though small, of finding the light quarks in the A;, baryon 
with longitudinal momenta of order Maj,. Therefore, reducible corrections on the A^ baryon 
side are dominated by soft dynamics, but contain weak double logarithms with collinear 
ones suppressed. Similarly, three-particle reducible corrections on the Ac baryon side are 
absorbed into the Ac baryon wave function. In the fast recoil region collinear divergences 
become stronger, and double logarithms associated with the Ac baryon wave function are 
more important. The remaining part of radiative corrections, with all collinear and soft 
divergences subtracted, are characterized by a scale of order Ma,,, and absorbed into the 
hard h quark decay amplitudes. Irreducible corrections, with a gluon attaching a quark in 
the At baryon and a quark in the Ac baryon, are infrared finite in the large recoil region 



|I3|, and also absorbed into the hard decay amplitudes. 

The kinematic variables are defined as follows. The A^ baryon is assumed to be at rest 
with the momentum 

]9=(p+,p-,PT) = ^(1,1,0). (4) 

The valence quark momenta in the A^ baryon are parametrized as 

ki = (p+, xip~, kix) , k2 = (0, X2P~, k2T) , k^ = (0, x-sp' , k^r) , (5) 

where ki is associated with the b quark. The momentum fractions and the transverse 
momenta obey the conservation laws, 

a;i + X2 + X3 = 1 , kiT + ^2t + k^r = ■ (6) 

The Ac baryon momentum is chosen as p' = {p'^ ,p'^ , 0) with p'^ ^ p'^ at fast recoil. We 
define the velocity transfer p, 



^ p-p' Ml + Ml 

Ma^ being the Ac baryon mass. Using the on-shell condition p' = M|^, the plus and minus 
components of p' are written as 

P'^ = P+P'^ , P'' = P-P' , (8) 

with 



p+ = {p+^p'-l)r, p_ = (p-v/p2-l)r, (9) 

and r = Ma^/Ma^. The valence quark momenta in the A^, baryon are parametrized as 

where A;'^ is associated with the c quark. The primed variables obey similar relations to 
Eq. (D. 

According to factorization theorem, the hadronic matrix element is expressed as 



M^= / [dx\[dx'] [d^kT][d^k'j.]^A.^,p,y{k[,fi)m''''^'''^^{k[,h,p,M^^,fi) 



X'^Aial3j{ki,IJ,) , (11) 

with the notations 

[dx] = dxidx2dx-i6 1 — X!^m ' [d'^kx] = d'^kiTd'^k2Td'^k3TS^ X! ^^^ ] ■ (-'-^) 

[dx'] and [d'^k'j] associated with the A^ baryon are defined in a similar way. The hard ampli- 
tude Hn will be computed in Sec. 111. The dependence on the factorization (renormalization) 
scale fi will disappear after performing a RG analysis. 

The structure of the A^ baryon distribution amplitude \l/ Aba/37 is simplified under the 
assumptions that the spin and orbital degrees of freedom of the light quark system are de- 
coupled, and that the A^ baryon is in the ground state (s-wave). The distribution amplitude 
is then expressed as 

6 



^ 



Abo/37 



(^i,/^) 



/At 



n^:^^^'"''^'^''(o|^K(^i)«J(^2)rf;(o)]|A,(p)) 



[(|^ + MAj75C]^^[Afc(p)],<l>(A;i,/i) 



(13) 



8V2Nc 

where iVc = 3 is number of colors, b, u, and d are quark fields, a, b, and c are color indices, 
a, P, and 7 are spinor indices, /a^ is a normalization constant, C is the charge conjugation 
matrix, and $ is the Af, baryon wave function. Under similar assumptions, the Ac baryon 
distribution amplitude '^\^ai3'y is written as 

1 



■[(;^' + MAj75C]/5.[Ac(p')].n(A;:,/i), 



/Ac 



(14) 



8V2N,' 

where the normalization constant f\^ and the wave function 11 are associated with the Ac 
baryon. 

Because of the inclusion of parton transverse momenta, Sudakov resummation for a 
hadron wave function should be performed in the impact parameter b space with b conjugate 
to A;t ^j^. The result is § 



$(fc. , bi, fi) = exp 



-J2s{w,ki)-3 — 7g(a,(/i)) 

1=2 •'^ ^^ 



(15) 



where 7^ = —as/ii is the quark anomalous dimension, and the factorization scale w is chosen 

as 



w = mm 



1 1 1 
61 ' 62 ' ^3 



(16) 



with 63 = |bi — b2|. The explicit expression of the Sudakov exponent s is given by |]T7 
s{w,Q) 



<3 dp 

w p 



ln(^JAiasip)) + B{aM) 



(17) 



where the anomalous dimensions A to two loops and B to one loop are 

) g / c7i5 \ / o, \ 2 

-nj + -(3o In (^ 



A = Cf— + 



67 TT^ 10 

Y ~ y ^ 27' 

,27B-1\ 



3 TT 



(18) 



Cp = 4/3 being a color factor, n/ = 4 the flavor number, and 7^; the Euler constant. The 
one-loop running coupling constant. 



TT 



/3oln(/iVA' 



(19) 



qcd; 



with the coefficient (3q = (33 — 2nj)/12 and the QCD scale Aqcd, "will be substituted into 
Eq. (p!7|). The initial condition of the Sudakov evolution absorbs nonperturbative dynamics 
below the factorization scale w. 

Following the derivation in P,|l8|, we obtain the Sudakov resummation for the Ac baryon 
distribution amplitude. 



n(A;-+,6i,/i) = exp 



, 1 Jw U 



tt(x',] 



(20) 



1=1 •'^ /^ 

We have included the Sudakov exponent s associated with the c quark, which carries large 
longitudinal momentum in the fast recoil region. Notice the same transverse extents hi as 
those for the Af, baryon. This is the consequence of neglecting the transverse momenta which 
flow through the virtual quark lines in H^ JTB|. 
The RG analysis of if^ leads to 



H^{k[^, ki , hi, p. Ma,, /i) = exp 



-3^/ — 7g(a,(/i)) 



H^{x[, Xi, hi, p. Ma,, ti, ta) , (21) 



where the superscripts a', /5', ■ ■ ■, have been suppressed. Since large logarithms have been 
collected by the exponential, the initial condition iJ^ of the RG evolution on the right-hand 
side of the above expression can be computed reliably in perturbation theory. To simplify 
the formalism, we shall make the approximations Mf, ^ Ma, and Mc ~ Ma^, and neglect 
the transverse momentum dependence of the virtual quark propagators as mentioned before. 
The two arguments ti and ^2 of H^, which will be specified in the next section, imply that 
each running coupling constant a^ is evaluated at the mass scale of the corresponding hard 
gluon. Substituting Eqs. ([T3|)-([2T|) into Eq. (|11]), we derive the factorization formula for the 
semileptonic decay A^ — > AJu, where the p dependence has disappeared as stated before. 



For the A;, baryon wave function (f){xi,X2,X3), we adopt the model proposed in |jT9 



0(C,r7) = iVr72^(l-r7)(l-C)exp 



M2 mf 



(22) 



2/32(1 -r^) 2/52^C(l-C). 
with A^ being a normalization constant, (3 a shape parameter, mi the mass of hght degrees 
of freedom in the A^ baryon. The new variables ( and rj are defined by 

C = r — ' V = X2 + X3 . (23) 

In terms of ( and rj, the normalization of </){(, rj) is given by 

'dCvdv<P{C,v) = l, (24) 



which determines the constant A^, once the parameters (3 and mi are fixed. The above wave 
function with the factor 1]'^C{1 — ri){l — () = a;ia;2a:3 suppresses contributions from the end 
points of momentum fractions. The exponents proportional to M^/(l —1]) = M^/xi and to 
m,'{/[ri({l — ()] = m^ /x2 + m^f/x^ with Mf, ^ m,i indicate that has a maximum at large 
Xi and at small X2 and 3:3, and that the b quark momentum kf is roughly equal to M^. For 
0(x3,xi,X2) which will appear in the factorization formulas presented in Sec. Ill, the above 
expression is transformed into 

M2 rnf(l-rj + rjCy 



(j){C,r]) = Nr]\{l-r]){l-C)exp 



(25) 



2/3277(1-0 2/3^r]Cil-v) 
For convenience, we assume that the Ac wave function 7r(0, rj') possesses the same functional 
form and the same parameters [3 and m,i as of (f){(, rj), but with the b quark mass Mj, replaced 
by the c quark mass Mc- The wave function tt{x[,x'2,x'^) also has a maximum at large x'l, 
such that the c quark momentum k^ is roughly equal to M^. 

III. TRANSITION FORM FACTORS 

In this section we present the factorization formulas for the form factors /i and gi, which 
are associated with the spin structures Ac7^Ab and Ac7^75A{, in M^, respectively. Working 
out the contraction of '^A^a'is'yH^ '^ '' °''^'^'^Atai3-y in momentum space, we extract the hard 
part H. Employing a series of permutations of the valence quark kinematic variables as in 



0], the summation over the leading diagrams in Fig. 1 reduces to two terms for each form 
factor. The factorization formula for the form factors /i(p) and gi{p) are written as 
47r f^ 1"°° ''^'"^ 



flip) = -7^ / [dx'][dx] / hdh,hdh2 / defU^ 
Zi Jo Jo Jo 



X J2 Hjix[, Xi, bi, p, MA^,tji)J^j{x[, Xi, p) exp[-S(x., x^, w, p, Ma„ tji)] , (26) 

5'i(p) = 7^ / [c?a;'][rfx] / bidbib2db2 / dOf^Jk^ 
27 Jo Jo Jo 

2 

X ^ /fj(Xi, Xi, bi, p, MA,,tji)gj{x[, x„ p) exp[-5(x., x^, w, p, Ma,, t^O] , (27) 

where 9 is the angle between bi and b2. 

The functions J^j and ^j, which group together the products of the initial and final 
baryon wave functions, are, in terms of the notations. 



0123 = 0(a;i, X2, X3) , 7ri23 = Tt{x\,X2, X3) 



(28) 



given by 
^1 



2 2Jp^-l-l){l-x[) 



!'i237ri23 [(1 - x[ - p-)p+ + r2](l - X[)X2P+ 

+ (2(1 + r)p - 4r - l)x2 + (2(2p - 1) + (2p - 3)pi) X2x[] 



+ 



[(1 - x[ - p-)p+ + r2](l - xi)a;2P+ 

(pi + 3)(1 - xi)x[ + 2 (2(p - l){,/^ + p) - A x'. 



pi + 2rJp^ - 1 + 3 + 4r - rp (1 - xi) 



(1 - Xi)2x2P+ 
r((2-p)pi + l-2p)(l-xi)x^] 
r 



2r 2Jp2 - 1 - 1 + 2p (1 - xi) + 2 Jp' - 1 - 2 + p x^ 



(1 -xi)(l -x'i)x2p+ 



2 (^V^P' - 1 + 2 - p) + r(pi + 3)(1 - xi)(l - xO 



-^2 



^312 



+2r(^2(p-l)(^p2-l+p)-ljx2 

r 

vr3i2 [(l-x[j -p_)p+ + r2](l -X3)xip+ 



(29) 



2rpi(l - x'^) + Ar^{l + p)(l - Xi) 



2 \ „/ 



+2r^ (3 - ^p2 _ 1 j a;^ _ 2r f (p - l)^p2 _ i - p^ j 3,; 

-(l+pi)(r(p-l)xi+x'2)(l-x;5)] 

2r 
+ [(x'2 - p-)(l - xi)p+ + r2][l - (1 - xip+)(l - x^)] 

10 



X 



r ip+ \J p^ - 1 j (X1X2 - Pi{xi + x'2)) +pi {/Ir'^xi + x'2 + 2r\J p^ - 1 

[2rpi(l-X3)+4(l + p)(l-x;) 



[1 - (1 - X2p+){1 - x'^)](l - Xs)p+ 
-r(l + P2) (^a;2 + ^J p'^ - lx[ ) (1 - X3) 



(30) 



^1 



^123 



vri23 [(1 - a;i - P-)p+ + r2](l - a;i)x2P4- 



[(2p-3 + (2p-l)p2)a:2(l-a:;: 



+2i2-p-Jp^-l]x2-2{2p-l + 2Jp^-l]{l-x[ 



2(2ip-l)i./^P^l+p)-l^x'2 



[(1 - x[ - p-)p+ + r2](l - Xi)x2P+ 
+2r (^V^P' - 1 + 1) (1 - ^1) - (3P2 + 1)(1 - a:i)(l - x'l) 

(2p - 3 + (2p - 1)P2)4(1 - a:i) + 2 (2 - p - Jp2 _ 1 ) ^'^ 



(1 - xiyx'2pi 



2i2p-l + 2Jp^-l]{l-x^ 



-2r(2(p-l)(v'p2-l + p)-l)x2 



(1 -a;i)(l -a;i)a;2P+ 
-2(v/p2 - 1 + 2 - p)(l - x[) + r(3p2 + 1)(1 - xi)(l - x'l) 



(31) 



rn , X , 2in ^' -4r2(l + p) + rp2(4 - p - p2)xi 

^31271-312 [(1 - 2:3 - p-)p+ + r^\ (1 - X3)a;'ip+ l 

+r(p - l)(p2 - 1)4x1 + 2r(l - 4) + 2r ("(p - l)(y'p2 - 1 + p) + l) 4 

-(p2 + l)4(l-4)] 

2r 
+ [{x'2 - p-)(l - xi)p+ + r2][l - (1 - xip+)(l - 4)] 

+r2xi + x'2 - 2r [p + y^p2 - 1] (1 - a;i)(l - 4) 



rK/pi-l-2-p 



[1 - (1 - a;2p+)(l - xi)\{l - x^)p+ 
-2 {jp^ - 1 + 2p - 1) 4 + 2r ((p - l)(v/p^-l + p) + 1) X2 

-r(p2 + l)(ra;2 + (p - 1)4) (1 " ^3)] , 



(32) 



with pi = y^(p + l)/(p - 1) and p2 = 1/pi. 
The hard parts are given by 



Hi = asitu)asiti2)Ko J {I - xi)(l - 4)p+^A6^i Kq J X2x'2P+Mi^Jj 



(33) 



11 



H2 = as{t2l)as{t22)KQ[^^^JxlX'lp+M^J)l\KQ[Jx2X'2P+M^J)2\ , (34) 

with Kq being the modified Bessel function of order zero. The complete Sudakov exponent 
S is written as 



S{x\, Xi, w, p, Ma^, tji) = Sd{x[, Xi, w, p, MaJ + Ss{w, tji) , (35) 



with 



3 3 

1=2 1=1 

Ss = 3 T' ^%{asm + 3 f' ^7.(«s(/i)) . (37) 

Jw p Jw p 

The hard scales tji are chosen as 



til = max 
^21 = max 



:1-Xi)(1-x'i)p+Ma„1/6i 



^Jxix[p+M/^^,l/bi 
t2i = t22 = max Jx2X2P+Ma^, 1/62 



(38) 



which are always greater than w. It is possible that the hard scales tji are small and the 
running coupling constants become large as bi are close to 1/Aqcd- These nonperturba- 
tive enhancements are, however, suppressed by the Sudakov exponential exp(— 5*^), which 
decreases quickly in the large 6, region and vanishes as hi > 1/Aqcd- The exponential 
exp(— ^d) approaches unity; that is, there is no Sudakov suppression from the all-order 
summation of infrared logarithmic corrections at small 6j. In these short-distance regions 



higher-order corrections are regarded as being hard and should be absorbed into H pO 
Another exponential exp{—Ss), as a consequence of single-logarithm summation, describes 
the RG evolution from the factorization scale w to the hard scales tji. 

For the case with massless leptons, it is easy to show that the differential decay rate in 
the rest frame of the A;, baryon is given by 



dT Ml 



3 



Afc' ^2 ITA |2 /„2 1 n f |2/„ iNro I Q„2 



G^lKfelVp' - l{l/ir(P - 1)[3 + 3r^ - 2(2p - l)r] 



dp 24n 

+ \9i\'{p + 1)[3 + 3r2 - 2(2p + l)r]} , (39) 

12 



where only the contributions from the form factors /i and gi are considered. It is straight- 
forward to obtain the total decay rate 

from Eq. ( PPJ ) and thus the branching ratio B{Ai, -^ AJu), if the form factors /i(p) and 
gi{p) in the whole range of p are known. 

IV. RESULTS 

In order to reduce the number of unknown parameters, we make an approximation. 
Consider the baryonic decay constant f\ defined, in heavy quark effective theory, by 

(0|?|Aq) = /aAq , (41) 



in terms of the A-baryonic current |21,22 



f = e'"'%u^C^,d'')K , (42) 

where Aq is the heavy baryon spinor, h^ the heavy quark field, and a, b, c denote the color 
indices. We contract a Dirac tensor (C75)/3^ with a heavy A-baryon distribution amplitude 
such as "^A^a/s-y in Eq. ([T3|) and integrate out the valence quark momenta /cj. Compared with 
Eq. (PI), we extract the baryonic decay constant 

/a = /acjMa^ . (43) 

It implies that in the heavy quark limit the normalization constants /a^ and /a^ are related 
by 

/a,Ma, = /a Ma, . (44) 

Therefore, f\^ associated with the Ac baryon will not be treated as a free parameter in the 
numerical analysis below. 
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We are now ready to compute the form factors /i(p) and gi{p) from Eqs. (^) and 
(p7|), adopting the CKM matrix element Vcb = 0.04, the masses M/^^ = 5.624 GeV and 
M\^ = 2.285 GeV, and the QCD scale Aqcd = 0.2 GeV. We examine the self-consistency 
of our calculation by considering the percentage of the full contribution to the form factor 
/i, that arises from the short- distance region with all as{tji)/iT < 0.5. The percentages for 
different (3 with rrii fixed at 0.3 GeV are listed in Table I. It is observed that the perturbative 
contributions become dominant gradually as p and P increase: a larger p corresponds to 
larger momentum transfer involved in decay processes, and a larger f3 corresponds to heavy 
baryon wave functions which are less sharp at the high ends of the momentum fractions Xi 
and x[. We conclude that the PQCD analysis of the transition form factors is self-consistent 
for P > 1.0 GeV and p > 1.2, viewing the perturbative percentage of about 80%. Compared 
to the corresponding meson decay B -^ Dlv 0, a perturbative expansion is less reliable in 
the baryon case, because partons in a baryon are softer, such that Sudakov suppression is 
weaker. 

To obtain the total decay rate, we need the information of /i and gi in the whole range 
of p. Since the perturbative analysis is reliable only in the fast recoil region, we extrapolate 



the PQCD predictions at large p to small p. Hinted by [^, we propose the following 
parametrization for the form factors: 

/i(p) = ^, 9i{p) = ^. (45) 

where the constants Cf and c^, and the powers af and Ug are determined by the PQCD 
results at large p. The constants Cf and c^, equal to the values of the form factors at zero 
recoil (p = 1), should be close to unity according to heavy quark symmetry. We fit Eq. (^Sj) 
to the PQCD results in the range with p > 1.3 for j3 = 1.0, where perturbative contribution 
has exceeded 80%. The powers aj = 5.18 and ag = 5.14, close to af ~ 4.6 at large p from 



the method of wave function overlap integrals [Q, are obtained. These values are larger 
than 1.8 extracted from the transition form factors associated with the corresponding meson 
decay B -^ Dlv [§]. This is expected, because perturbative baryon decays involve more hard 
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gluon exchanges. 

On the experimental side, there exist only the data of the semileptonic branching ratio 
-B(Ab -^ XW) ~ 10% 1^, where the final-state particles X are dominated by the charm 
baryons. The data of the B meson semileptonic decays show B{B -^ D*W) ~ 3B{B — »• 
Dlu), indicating that each of the three polarization states of the D* meson contributes the 
same amount of branching ratio as the D meson does. It is possible that this observation 
applies to dominant modes in the Af, -^ XW decays with the excited charm baryons Ac(2593) 
of spin J = 1/2 and Ac(2625) of J = 3/2. That is, the branching ratio B{Ab — ^ Ac/z/) is about 
1/4 of B{Ah — > XW), i.e., about 2 ~ 3%. This estimation is consistent with the experimental 
upper bound of the branching ratio from the data B{Ah -^ AJu + X) = (8.27 ±3.38)% p5[ . 



We substitute Eq. (^) for the form factors /i and gi into the decay rate T in Eq. (^OD, 
and adjust the normalization constant f\^ such that our predictions for the branching ratio 
are located in the range of 2 ~ 3%. The Ac baryon normalization constant /a^ changes 
according to Eq. (||). We adopt the A^ baryon lifetime r = (1.24 ±0.08) x 10"^^ s ^. The 
value of /aj determines the parameters Cf and Cg. It is then found that /a^ = 2.71 x 10~^ 
GeV^, corresponding to 

. / ^ 1-32 , , -1.19 

gives the branching ratio 2%, and /a^ = 3.0 x 10"'^ GeV^, corresponding to 

, , ^ 1-62 , , -1.46 ,^^, 

/i(/^) = -JT?, ' Sx{p) = ^JT ' (47) 

gives the branching ratio 3%. Since the values of the form factors at zero recoil should be 
close to unity as stated above, we prefer Eq. (|46|) with /i(l) = 1.32 and gi{l) = —1.19, which 
are also consistent with the conclusion in ||2^. The corresponding normalization constant 
/a, = 2.71 X 10-3 GeV^, of the same order as fp = (5.2±0.3) x lO^^ GeV^ for the proton @, 



is reasonable. The PQCD predictions and the corresponding extrapolations are displayed in 
Fig. 2, which deviate from each other at small p. If applying the PQCD formalism to the 
zero recoil region, we shall obtain divergent form factors as shown in Fig. 2, which imply 
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the failure of PQCD. Note that our results of the form factors exhibit slopes larger than the 



dipole behavior assumed in ||23|| . 

We then examine the sensitivity of our predictions for the branching ratio -B(Afe — »• 
Ac/z/) to the variation of the parameter /3. Choosing /5 = 2.0 GeV and (3 = 4.0 GeV, and 
normalizing the corresponding form factors in the way that they have similar values to those 
for /3 = 1.0 GeV in Eq. (|46|) , we obtain the form factors 

1.34 , . -1.17 



flip) 



5.04 



and 



flip) 



P 



1.34 



9iip) 



,4.92 



9iip) 



P 



-1.18 

o4.79 



(48) 



(49) 



p^.^-* p^ 

respectively. Equations (|48|) and (|^) lead to increases of the branching ratio by 4% and 
8%, respectively. That is, our predictions for the branching ratio are not sensitive to the 
choice of baryon wave functions. This observation is attributed to the fact that the PQCD 
results of the transition form factors at large recoil are insensitive to the variation of baryon 
wave functions. 

We present in Fig. 3 the differential decay rate dV /dp derived from the form factors in 
Eq. (^61) , which can be compared with experimental data in the future. The A;, and Ac 
baryon wave functions determined in this work are given by 



m. 



0(C,r7) = 6.67x10^2 ^2^(1 _^)(1_^) 

r M,^ . 

^""""^L 2(1.0 GeV)2(l-r7) 2(1.0 GeV)2r/C(l - C) 
7r(C,r/) = 6.94xl0%2^(l-r7)(l-C) 



X exp 



Ml 



mf 



(50) 



(51) 



2(1.0 GeV)2(l-r/) 2(1.0 GeV)2r/C(l - C). 

At last, we compare our predictions with those derived from other approaches in the 

literature. The A^ — > A^ transition form factors have been evaluated by means of overlap 

integrals of infinite-momentum-frame (IMF) wave functions, nonrelativistic and relativistic 
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quark models, and QCD sum rules. For a review, refer to [^. Basically, they are nonpertur- 
bative methods without involving hard gluons. QCD dynamics is completely parametrized 
into IMF wave functions in the overlap- integral approach [^,^], and into baryon-three- 
quark vertex form factors in the relativistic quark model ||29|. Information of the above 
bound-state quantities can be obtained by solving Bethe-Salpeter equations PD[- Most of 



the analyses, including QCD sum rules |2^j31|j3^ , led to the branching ratios about or below 



6%. The prediction -B(A{, -^ AJiy) ~ 9% in p8| is a bit higher compared to the data of 



B{Ab -^ AcW + X). Our result is close to (3.4 ± 0.6)% derived in fp. 

V. CONCLUSION 

In this paper we have developed PQCD factorization theorem for the semileptonic heavy 
baryon decay Af, -^ AcW, whose form factors are expressed as the convolutions of hard b 
quark decay amplitudes with universal A^ and A^ baryon wave functions. It is observed that 
the PQCD formalism with Sudakov suppression in the long-distance region is applicable to 
Af, -^ Ac decays for the velocity transfer greater than 1.2. This observation indicates that 
PQCD is an appropriate approach to analyses of two-body exclusive nonleptonic A^ baryon 
decays. Requiring that the normalizations of the form factors at zero recoil are consistent 
with heavy quark symmetry, we have predicted the branching ratio B{Ai, -^ AJi?) ~ 2%. 
We have also determined the A^ and Ac baryon wave functions shown in Eqs. (^Oj) and (pl|), 
respectively. These wave functions, because of their universality, will be employed to study 
nonleptonic A;, baryon decays in the future. 

This work was supported by the National Science Council of the Republic of China under 
Grant Nos. NSC-88-2112-M-001-041 and NSC-88-2112-M-006-013. 
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FIGURE CAPTIONS 

FIG. 1 Lowest order diagrams for the A5 -^ AJi? decay. 

FIG. 2 Dependence of /i and \gi\ on p for j3 = 1.0 and mi = 0.3 obtained from PQCD 
(solid lines) and from the extrapolation in Eq. (|46|) (dashed lines). The upper (lower) set of 
curves represents the form factor /i (1(711). 

FIG. 3 Dependence of dV /dp on p obtained from Eq. (^) in units of lO^^'^ GeV. 
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Table I. Percentages of perturbative contributions for various f3 and p. 



percentage 


p= 1.2 


p = 1.3 


p = 1.4 


f3 = 1.0 GeV 


77.7% 


83.6% 


85.2% 


(3 = 2.0 GeV 


79.3% 


83.0% 


85.7% 


(3 = 4.0 GeV 


82.3% 


84.7% 


86.3% 
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